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Abstract 

The coefRcients of dispersion, Dij, and the dispersivity, Qijki, appear in the expression for 
the flux of a solute in saturated flow through porous media. We present a detailed analysis of 
these tensors in an axially symmetric porous medium and show that in such a medium, the 
dispersivity is governed by 6 independent moduli. We present also the constraints that have 
to be satisfied by these moduli. We also show that at least two independent experiments are 
required in order to obtain the values of these coefficients for any three-dimensional porous 
medium domain. 

1 Introduction 

The coefficients of dispersion and dispersivity appear in the expression for the flux, J, of a solute 
in saturated flow through porous media. Here, a porous medium domain is one that is made up of 
two parts: a solid matrix and a void-space completely occupied by one or more fluid phases. The 
void-space is interconnected, thus allowing fluid flow through it [3J. Both parts are distributed over 
the entire porous medium domain. We study the transport of a solute dissolved in a fluid phase. 
At every point within the fluid, a solute concentration, 7, and a fluid velocity, V, can be identified. 
A mathematical model describing the solute's transport within the void space requires information 
on the fluid-solid interface that bounds the fluid occupied domain. However, in natural materials, 
such as soil and fractured rock, this information is not available. Therefore, the problem of solute 
transport is usually described in terms of averaged velocity, V, and averaged solute concentration, 
7. These averaged values, taken over some representative elementary volume, are assigned to every 
point in the porous medium domain. 



With the total (advective) solute flux at a point inside the void-space defined by 7V, and at a 
point in the porous medium domain, regarded as a continuum, defined by 7V, we have j^. 



7V = 7V + 7'V', where 7' = 7 - 7, V' = V - V, 

are deviations from the average values, with 7' = V = 0. The above equahty states that the total 
flux of the solute at a point in the porous medium domain, is made up of an advective flux (= 
product of average concentration and average velocity), and another flux, called dispersive flux. 



equal to the average of the product of the two deviations, J = 7'V'. 

The above discussion can be extended to any extensive quantity, with 7 representing its intensive 
quantity. Thus, the phenomenon of dispersion occurs whenever an extensive quantity (mass of a 
solute, in the example considered above) is transported in a fluid phase that flows in the void space 
of a porous medium. It is a consequence of using an average velocity and average concentration to 
describe the advective flux of the solute in the void space. The total flux is the sum of the advective 
flux and the dispersive one. In this article, we use the mass of dissolved chemical species as the 
considered extensive quantity. For short, we denote the average values: 7 = c, V = U. 

The coefficient of dispersion, Dij, is a symmetric tensor of 2nd rank, which appears in the 

Fickian-type expression for the dispersive flux of a solute (mass of solute per unit area of void 

space in a cross-section of the porous medium) [3] , 

3 

Ji = - ^ DijVjC, Dij = Dji, i, j = 1,2,3 (1) 

where Jj denotes the ith. component of the solute's 3-dim flux vector, J, and c is the solute's average 

concentration. Equation ([1]) is valid for the general case of any anisotropic porous medium, with 

isotropic media as a special case. 

In contrast to similar linear flux laws, e.g.. Ohm's and Fourier's laws, with corresponding 

tensorial coefficients of electroconductivity, aij, and of thermoconductivity, Xij, respectively, which 

depend only on properties of the considered medium, here, in the case of the dispersive flux in 

saturated flow through a porous medium, the dispersion tensor Dij depends also on the fluid's 

average velocity U in the void space 

^3 3 
= -[jYl = ^ Ui, U^ = U^, (2) 

k,l=l k=l 

where aijki is a property, called dispersivity, of the porous medium only, and Ui is the ith component 
of the fluid's average velocity vector U. The dispersivity, is a 4th rank tensor, with intrinsic 
symmetry [^^]^ in Jahn notations [5] 

[V'^] : Qijki = ttijik , aijki = ajiki- (3) 



Additional simplifications in dealing with the dispersivity tensor can be obtained by making 
use of certain auxiliary algebraic relations among its components. These relations appear under 
the action of the crystallographic point symmetry group G C 0(3) of the porous medium that 
makes some (not all) of the 36 elements of aijki survive, or at least independent of the remaining 
ones. According to [9], the decrease in the number N{Q) of a^fci-components is significant when 
we increase the symmetry, i.e., 

N{Ci) = 36 ^ ... ^ N{V2h) = 12 ^ N{Vih) = 7 ^ N{V^h) = 6 ^ A^(C(3)) = 2. (4) 

In thermodynamics, the rate of entropy production, o", is related to the thermodynamic driving 
force, X, and to the thermodynamic flux, Y. According to [4J, two vectors X and Y are considered 
conjugated if they obey the relationship 

3 

a = Y,YiX,. (5) 

2 = 1 

In the case considered here, the dispersive flux of the solute, J, is driven by Vc, which acts as a 
'driving force.' In other words, X = gVc and Y = J, where g denotes a (dimensional) scalar 
parameter, which depends on the considered transport phenomenon. In this case, the rate of 
entropy production is expressed by 

3 

By requiring that the quadratic form Q be positive definite, we derive the following constraints 
on the Djj-matrix, 

f Dii Di2 1 f Du Di3 ] f D-?-? D23 1 

Da > 0, det ^ } > 0, det <^ } > 0, det ^ " V > 0, det Aj > 0. 

[ D21 D22 J [ I?31 I?33 J [ D23 D33 J 

The above constraints hold irrespective of the fluid's averaged velocity U. Furthermore, by ([2]), 
they provide another set of constraints imposed on elements of the dispersivity tensor, aijki- Such 
straightforward way to establish these constraints is not effective, since it requires a lot of routine 
algebra and does not prevent repetitive inequalities. In section HI we shall present a much more 
effective method for deriving the sought constraints. 

In the meantime, let us estimate the total number of such inequalities. Making use of the 
2-dim matrix representation of the aijki-tensor [9\, one can show that the number M{Q) of such 
constraints decreases also with growth of the symmetry Q of the porous medium, 

M(Ci) = 63 ^ . . . ^ M{V2h) = 10 ^ M(p4h) = 7 - M{V^h) = 6 ^ M(0(3)) = 2. (7) 
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2 The Tensor aijki in a Porous Medium with Uniaxial Symmetry 

In this section, we study the 3-dim tensor aijki whose components depend only on porous medium 
properties, actually, for saturated flow, only on the geometry of the medium's void space. We shall 
focus on porous media with uniaxial symmetry. Our interest in this kind of porous medium is 
motivated by two reasons. First, in saturated flow in ground water geological formations, called 
aquifers, the porous medium comprising the latter is often axially symmetric [8j, with a vertical 
axis, e.g., in a layered aquifer. On the other hand, four decades of continuing discussions (see |7J 
and references therein) on explicit expressions of the tensors aijki and Dij in porous media with 
uniaxial symmetry have not provided any rigorous analysis of these entities. 

Instead of a 2-dim matrix representation of the Ojjfci-tensor, presented in [9], we suggest an- 
other representation, which originates in the theory of uniaxial nematic liquid crystals (LC). The 
underlying idea is that when we compare the ajjfc/-tensor with the 4th rank viscosity tensor, rjijki, 
of a uniaxial LC, its intrinsic symmetry is [^^]^ , i-e., slightly 'stronger' than ([3]), with, 

[V^] : rjijki = r]ijik, ijij^ = rjjiki, rjij^i = rjkuj. (8) 

According to the known representation of rjijki is 

'nijki = m^ij^ki + m (SikSji + SuSjk) + -nu {ei^jSui + e^eiSij) 

+i]5 {eiCkSji + ejCkSii + eiei5jk + ejeiSik) + rj^eiejekei, (9) 

in which the e^'s denote the three components of the unit vector e, along the symmetry axis of the 
medium. 

From ([9]), it follows that there exists only one term, {eiej5ki + ekei6ij), which reduces its intrinsic 
symmetry: [^^]^ [y^] ^, after decomposition of this term into two separate terms: ejCjJfc/ and 
Gk^l^ij- This explains why, in accordance with ([7]), the tensor a^jki has six independent moduli and 
how to construct its most generic form: 

aijki = ai 6ij6ki + i^ik^ji + hi^jk) + ^3 eiCjdki + 04 ekei5ij + 

— {eiCkdji + ejCkSii + CiCidjk + ejCiSik) + eiejCkei. (10) 

Note that by making ^ in (jlOp . we obtain the dispersivity tensor for an isotropic porous 
medium, 

a-ijki = aiSijSki + Y {SikSji + SuSjk) ■ (11) 

The axisymmetric tensor (|1U|) has to be supplemented by a set of constraints that results from 
the thermodynamic requirement imposed on the rate of entropy production, a, defined in ([6|). 



3 The Tensor Dij in a Porous Medium with Uniaxial Symmetry 



In this section, we consider the 3-dim tensor Dij and certain interesting related phenomena. Based 
on ([2]) and (llOh . we rewrite this tensor as 



UDij = aiU'^Sij + a2UiUj + asU^ejej + 04(6, U)^(5ij + 05(6, U)(C/iej + eiUj) + 05(6, U)^ejej, (12) 

where (e, U) denotes a scalar product of the two vectors: e and U. 

It is interesting to note that Batchelor [1], in his work on axisymmetric turbulence, derived a 
general expression for Dij that is based on five tensorial terms: CiCj, 6ij, UiUj, UiCj, UjCi, and four 
scalar functions, Cj(U,e), 



UDij = Ci{lJ,e)eiej + C2{V,e)6i, + Cs{V,e)UiUj + ^(U, e)(C/ie, + UjCi). 



(13) 



Keeping in mind that there exist only two quadratic-in-C/ scalar invariants, [/^ and (U, e)^, and 
only one linear- in- {/ scalar invariant, (U,e), we can verify the equivalence of the two expressions, 
([12]) and ([13]), by noting that 



Ci(U,e)=a3U^ + ae{U,ef, C2(U , e) = a^U'' + a^{\J , ef , C3(U,e) = a2, C4(U, e) = a5(U, e). 



We introduce the two 6-dim vectors, Dg and ag. 



De = {Du, D22, D12, D23, D31} , ag = {01,02,03,04,05,06}, 



(14) 



and denote their transposed counterpartners by Dg and ag . Next, we rewrite (|12p in the form 

1 



D 



g -- Ai(U,e)-a^, 



(15) 



where the 6x6 matrix Ai (U, e) is defined as 



Ai(U,e) = <^ 



2(e,U)eiC/i 
2(e,U)e2C/2 
2(e,U)e3t/3 



(e,U)2ef 
(e,U)2ei 
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(e,U)2ei 



(16) 



U2 Uf efu2 (e,U)2 

U2 Ul eiU2 (e,U)2 

U2 Ui eiU2 (e,U)2 

U1U2 eieaU^ (e, U)([/ie2 + eiC/2) (e,U)2eie2 

U2U3 eaesU^ (e, U) ([7263 + 62^/3) (e,U)2e2e3 

[ U3U1 eseiV (e,U)([/3ei + e3?7i) (e,U)2e3ei J 

Equation (jlSh opens the way to the derivation of the 6 dispersivity Oj-moduli from values of 
Dij. The latter can be obtained from experimental measurements, assuming that the symmetry 
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axis e is known, and so is the velocity U, in an appropriate reference frame, xi,X2,X3. However, a 
straightforward calculation shows that the matrix Ai (U, e) is degenerated, 

det Ai (U,e) = 0, rank Ai (U, e) = 4 . (17) 

Equality (llTh leads to two most important consequences and raises a question: 

1. At most, four independent components Dij can be obtained from a single experiment. The 
remaining two -Dj^'s are linearly representable by the first four. 

2. One cannot get all six dispersivity aj-moduli from a single experiment. 

3. How many experiments do we have to perform in order to obtain all six Oj-values? 

In what follows, we shall demonstrate that two different experiments suffice for obtaining all six 
Oj-values. We consider two cases (= experiments): 

1. Ui = 0, U2 = Us = u and ei = 1, 62 = 63 = 0. Then, 

Du = 2(ai + a3)u, D22 = D33 = (2ai + a2)u, D23 = a2U, Du = D-n = 0. (18) 
The above relationships provide three moduli: ai, 02 and 03. 

2. Ui = U2 = u, U3 = Q and ei = 1, 62 = 63 = 0. Then, 

Dii = (2ai + 02 + 203 + a4 + 205 + ae)^, -D12 = (02 + a^)u., (19) 
D22 = (2ai + 02 + a4)n, -D33 = (2ai + 04)^, -D23 = -D31 = . 

The above relations provide three additional moduli: 04, 05 and a^. 

Thus, by two experiments, with different setups of U, and a known e, we have obtained the 
entire set of 6 dispersivity moduli in a 3-dim saturated porous medium with uniaxial symmetry. 

3.1 The Tensor Dij in an Isotropic Porous Medium 

For an isotropic porous medium, instead of ()12p . we have 

UDij = aiU^Sij + a2UiUj. (20) 

This means that only two moduli: ai and 02 are required for a full description of the Dij- 
components. The expression for Ai (U, e)-matrix, which is defined in (jl6p . is even more degener- 
ated, with rankAi (U, e) = 2. This indicates a strong linear dependence among the Djj-values: 
there exist four linear relations among the six Djj-elements. 
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Let us consider an isotropic porous medium, with the setup: Ui = u, U2 = U3 = 0. Then, 




(21) 



In ground water hydrology [3], the modulus oi is called transversal dispersivity (denoted ax) 
and oi + a2 is called longitudinal dispersivity (denoted ai)- 



4 Thermodynamic constraints 

In this section we determine the algebraic constraints that have to be satisfied by the six dispersivity 
Oj-moduli for an axisymmetric porous medium. We substitute (fT2]) into ([6]) and present a in its 
invariant form, irrespective of the reference frame. 



= aiU'(Vc)" + a2 (U,Vc)' + a3U'(e,Vc)' + a4 (e,U)'(Vc) 
+2 as (e, U) (U, Vc) (Vc, e) + ae (e, U)^ (e, Wcf . 



(22) 



We then represent (I22p as a quadratic form. This is a standard way to write the necessary and 
sufficient conditions for a to be positive definite. For this purpose, we choose the reference frame 
such that ei = 62 = and 63 = 1. By substituting these into (j22]) . we get 



i=l i=l 



i=l 



4=1 



+2a5f/3V3C^C/iViC + a6C/|(V3c)2 . 



(23) 



i=l 



Expression (|23|) is a quadratic form in a 9-dim space, with orthogonal basis UiVjc, i,j = 1,2,3. 
We decompose it in the form 



a 



72 = "1 



C/i (V2c)2 + Ul [Vicf] + (ai + 03) (V3C)' [C/f + + (ai + a4)Ui \{Vicf + (V2C) 



+ (01+02) Uf{ViCy + UiiV2Cy +(0i + 02 + 03 + a4 + 2a5 + 06)f/|(V3c)^ 

+ 2a2i7iC/2VicV2C + 2(02 + 05)^/3 V3c(;7iVic + C/2V2c) , 
and require the positive definiteness of the last expression. Then we obtain three inequalities 



(24) 



oi > , oi + 03 > , oi + 04 > , 



(25) 
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and requirement of positive definiteness of the 3x3 matrix T, 



ai + 02 02 02 + 05 

r = < 02 oi + 02 02 + 05 > • 

02 + 05 02 + 05 Oi + 02 + 03 + 04 + 2O5 + 06 

The last claim leads to three additional algebraically independent restrictions, 

oi + 02 > , oi + 02 + 03 + 04 + 205 + 06 > , 

of + oi(3o2 + 03 + 04 + 205 + og) + 202(03 + 04 + oe) > 2og 



(26) 



(27) 



Thus, we have arrived at 6 constraints imposed on the 6 dispersivity Oj-moduli; this is in full 
agreement with ([7]) for uniaxial symmetry group Q = Pco/i- Note that the first two inequalities in 
([25]) and (p7|) correspond to isotropic porous media. 



5 The Tensor Dij in a Planar Porous Medium with 
Dihedral Symmetry 

In this section we consider the dispersion phenomenon in a 2-dim porous medium domain, which 
is obtained by intersecting a 3-dim uniaxial porous medium with a plane that is parallel to the 
uniaxial axis. In the resulting 2-dim porous medium domain, there are two mirror lines, which are 
perpendicular to each other. By these symmetry elements, it is easy to recognize the dihedral point 
symmetry group 1)2, which is one of ten 2-dim point symmetry subgroups of 0(2). 

Prom the mathematical standpoint, the dispersion tensor, Dij, can be represented by a 2 x 2 
symmetric matrix, with three components: D3 = {Dn, D22, D12}, instead of the 6x6 matrix 
A2 (U, e), as in (fTC|) . By doing so, we obtain: 



D 



1 



^ A2 iV,e)-4, 



(28) 



A2(U,e) = <^ 



(e,U)2e? 
(e,U)2ei 



(29) 



where the 3x6 matrix, A2 (U,e), is defined as 

U2 Uf efU2 (e,U)2 2(e,U)eiC/i 
U2 [7| e^U^ (e,U)2 2(e,U)e2C/2 \c, v./ 02 

U1U2 eie2U2 (e,U)([/ie2 + eiC/2) (e,U)2eie2 

Equation (j28p does not enable the determination of the six dispersivity Oj-moduli, from 3 values 
of Dij, which can be obtained from a single experiment. For comparison, in section O we have 
shown that by two experiments, with different setups of U, and a known e, we have obtained the 
entire set of 6 dispersivity moduli in a 3-dim saturated porous medium with uniaxial symmetry. 



In this conjunction, the fohowing question arises: how many experiments do we have to perform 
in a 2-dim saturated porous medium domain with dihedral symmetry, in order to obtain the entire 
set of 6 dispersivity moduh? 

It is easy to show that two experiments are not enough. Indeed, let us consider two different 
flow setups, with fluid averaged velocities W and U, and the same orientation of e, 

iDl^ = i A2 (W, e).^J , 2Dl^ = i A2 (U, e) • a^^. (30) 

In (j30p the front subscripts 1 and 2 in the notations 1D3 and 1 indicate data obtained from the 
1st and 2nd experiments, with velocities W and U, respectively. We now introduce a new 6-dim 
vector, D6,a, 

De.a = {iDn, 1D22, iDu, 2Du, 2D22, 2D12}, (31) 

composed of two 3-dim vectors 1D3 = {i-Dii, 1-D22, iDu} and 2D3 = {2D11, 2D22, 2-^12}, and 
rewrite the two equations in (j30p . in the form 

D^, = l A3(W,U,e).a^, where (32) 



A3 (W,U,e) = <^ 



W2 






(e,W)2 


2(e,W)eiVFi 


(e,W)2e? 




wi 




(e,W)2 


2(e,W)e2W2 


(e,W)2ei 





WIW2 


6162 





{e,W)iWie2 + eiW2) 


(e, W)2eie2 




ul 




(e,U)2 


2(e,U)ei?7i 


(e,U)2e? 








(e,U)2 


2(e,U)e2f/2 


(e,U)2ei 





UIU2 


eie2U^ 





(e,U)([/ie2 + ei?72) 


(e,U)2eie2 



A straightforward calculation shows that 

det A3 (W, U, e) = 0, rank A3 (W, U, e) = 5. (33) 

Thus, among the six elements of the vector Dg^ai given in (|3ip . there are only five independent 
ones, and equation (j32p does not suffice for providing the six dispersivity aj-moduli. 

We shall now show that three experiments are sufficient for providing the entire set of a^'s. 
Consider 3 different flow setups, with velocities: W, U and T and with the same orientation of e, 

iD^ = i A2 (W, e).4 , 2D^ = i A2 (U, e) • a^, , 3^1 = ^ ^2 (T, e) • a^, (34) 

where 1D3 and 2D3 were defined earlir and 3D3 = {^Dn, 3D22, 3-D12}. In contrast to ()3ip . we 
introduce another 6-dim vector, Dq i,, 

D6,6 = {iDn, iDi2, 2Dn, 2D12, 3^11, 3^12}, (35) 
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composed of two elements, {Du and iDi2, of three 3-dim vectors jDa, i = 1,2,3. We rewrite the 
three equations ()34p in the form, 



Di;, = i A4(W,U,T,e)-a^, where 



(36) 



A4(W,U,T,e) = <^ 



W2 








W1W2 eieaW^ 
U1U2 eieaU^ 



1T2 6162!^ 



(e,W)2 


(e,U)2 


(e,T)2 




2(e,W)eiVFi (e,W)2ef 

(e,W)(VFie2 + eiW2) (e,W)2eie2 

2(e,U)eiC/i (e,U)2ef 
(e,U)(C/ie2 + eiC/2) 

2(e,T)eiTi 
(e,T)(Tie2 + 61X2) 



(e,U)2eie2 
(e,T)2eie2 



Note that the three 2-dim vectors: W, U and T are always coplanar. However, if at least two 
of them would be collinear, then ([36]) cannot be solved uniquely, i.e. det A4 (A:U, U, T, e) = 0, 
where A; is a real number. Otherwise, a straightforward calculation shows that if neither two of the 
three vectors: W, U, and T are collinear, then 



det A4 (W, U, T, e) / 0, rank A4 (W, U, T, e) = 6. 



(37) 



Thus, three special (not any) experiments suffice for providing all six dispersivity Oj-moduli in a 
2-dim saturated porous medium with dihedral symmetry. 



6 Concluding Remarks 

For the case of an axisymmetric porous medium, we have shown that six independent aj-moduli 
are needed in order to determine all components of the dispersivity tensor, and the constraints 
that these moduli have to satisfy. This information is required when determining the latter by 
experiments for specific porous media. We have also found the number of experiments that are 
required to determine the entire set of dispersivity moduli for 2-dim and 3-dim domains. This is 
important for experimental determination of the dispersivity. 
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